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Abstract 

For systems of N charged fermions (e.g. electrons) interacting with longitudinal op- 
tical quantized lattice vibrations of a polar crystal we derive upper and lower bounds on 
the minimal energy within the model of H. Frohlich. The only parameters of this model, 
after removing the ultraviolet cutoff, are the constants U > and a > measuring the 
electron-electron and the electron-phonon coupling strengths. They are constrained by the 
condition ^/2a < U, which follows from the dependence of U and a on electrical prop- 
erties of the crystal. We show that the large TV asymptotic behavior of the minimal energy 
Em changes at \/2a — U and that \/2a < U is necessary for thermodynamic stability: 
for \/2a > U the phonon-mediated electron-electron attraction overcomes the Coulomb 
repulsion and £jv behaves like —N 7 ^ 3 . 
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1 Introduction 

We study a system of N electrons in a polar (ionic) crystal, modelled by a Hamiltonian derived 
by H. Frohlich lfl2l . The model takes into account the electron-electron Coulomb repulsion, 
and a linear interaction of the electrons with the longitudinal optical phonons. The model is 
called the 'large polaron' model, since it assumes that a polaron (dressed electron) extends over 
a region which is large compared to the ion-ion spacing. In particular the underlying discrete 
(and infinite) crystal is replaced by a continuum. See lF71 fTTlfT9Tl . 

As is well-known, linear electron-phonon couplings induce an effective pair attraction be- 
tween electrons. This attraction competes with the electron-electron repulsion and may cause 
a phase-transition as the electron-phonon interaction strength increases. This mechanism is 
behind the production of Cooper pairs in the BCS model of low temperature superconductiv- 
ity, and in high-T c superconductivity the role of many-polaron systems is being investigated 

The Frohlich Hamiltonian depends on two non-negative dimensionless quantities, U and 
a. The constant U is the electron-electron repulsion strength, and a is the Frohlich electron- 
phonon coupling constant. Physically relevant models must satisfy the constraint, cf. B1I281. 

y/2a < U. 

In this paper we prove upper and lower bounds on the minimal energy En of the N- 
electron Frohlich Hamiltonian for all N and all non-negative values of U, a. In the unphysical 
regime our results imply that En ~ —N 7 / 3 . In the physical regime we find 

that En > —CN 2 , thus establishing a sharp transition in the large iV-asymptotics of En at 
y/2a = U. This transition is due to the mediated attraction between electrons overcoming the 
repulsion at y/2a = U in the limit of large N. In fact, the quantity U — V2a appears in our 
analysis as an effective Coulomb coupling strength. We also demonstrate that En < —aN and 
En+m < En + Em in the physical regime. We do not know whether or not En is an extensive 
quantity, but if it is not extensive, then this must be due to electron-phonon correlations, cf. 
Proposition IA.3I 

We pause this discussion to introduce the mathematical model. The Frohlich Hamiltonian 
describing N electrons in a polar crystal reads 

N 

J2[-^A x , + ^(x e )]+H ph + UV c , (1.1) 
i=\ 

where the number operator 

H p h = / a* (k)a(k)dk, 
accounts for the kinetic energy of the phonons while the field operator 

$( s )= f -J—\ e ^a(k) + e- lk - x a*(k)]dk, 
Jrz co\k\ 

is responsible for the electron-phonon interaction. Here cq := 2 3 / 4 7r. Finally the electron- 
electron interaction is given by the sum of two-body Coulomb potentials 

Vc(xi,...,x N )= V] 1 _ r 

l<i<j<N 1 ■>' 



N-Polaron Systems 



3 



We work in units where the frequency of the longitudinal optical phonons, wlo> Planck's 
constant h, and the electron band mass are equal to one. 

Let T denote the symmetric Fock space over L 2 (R 3 ). The Hamiltonian (11.11 ) defines a 
symmetric quadratic form onH = A N L 2 (R 3 ) but, a priori, it is not well defined as a self- 
adjoint operator. For that one must first impose an ultraviolet cutoff on the electron-phonon 
interaction: Let A > 0, and define the cutoff Hamiltonian as 

N N 

e=i e=i 

where 

$ A (x)= / AnW k ' x a(k) +e- lk - x a*(k)]dk. 
J\k\<A co\k\ 

The operators i?j\r,A are self-adjoint on P^phjnP^^ A^), by the Kato-Rellich theorem, 
and it is well known, cf. J2j |6j [13] EJ [271, that i?jv,A converges, as A — > oo, in the norm- 
resolvent sense to a self-adjoint operator, which we denote by Hjy. This implies that 

E N = lim E NA (1.2) 

A^oo 

if Ejy,A ■= inf cr(i/jv,A) an d En := inf a(H]y). 

The main goal of this paper is to investigate the large N behavior of the minimal energy 
En as a function of a and U. Our first result is an upper bound in the regime 

Theorem 1.1. There is a constant C such that for all N and for \[2a > U > 

E N < (V2a - Ufm [£ PTF + CN~vf] . 
Here Eptf < is given by (11.31) below. 

Theorem 11.11 is proved variationally by using Pekar's ansatz in terms of a product state, 
which is known to give the correct ground state energy for N = 1,2 in the large a limit 
(9l[25l[26l. Taking the expectation value in a state / <g) j] € A N L 2 (R 3 ) ® T and explicitly 
minimizing with respect to -q we arrive at a Hartree-Fock type energy which is then estimated 
by a Thomas-Fermi energy. This allows us to scale out all parameters and we are left with the 
bound in Theorem ll.il where 

-Eptf = inf £ PT fO), (1-3) 

p>0,J p(x)dx=l 

£ptf(p) := ^(6vr 2 )l f p(x)Ux - \\ ^IMdxdy. (1.4) 
Jr3 7 M e \x - y\ 

We note that in the error term in Theorem ll. li the exponent 1/17 can be replaced by any number 
less than 2/33 at the expense of a larger and divergent constant C. 

To show that the variational upper bound from Theorem ll. II has the right asymptotics in N 
and a, we provide the following lower bound: 

Theorem 1.2. There exists C > such that for all N and \[2a > U > 0, 

E N > -C G {V2a-U) 2 m -Ca 2 m~K (1.5) 
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This lower bound is obtained, essentially, by completing the square with respect to cre- 
ation and annihilation operators in the expression H p h + \/^Sj=i ^( x j)- The computation 
brings out an effective Coulomb interaction with coupling strength — y2a. Unfortunately, it 
also yields an infinite self-energy, which must be dealt with before completing the square. For 
that we use a commutator argument from [25], which is responsible for the error term in Theo- 
rem [L2] The resulting effective Hamiltonian with an attractive Coulomb potential is bounded 
below by the 'gravitational collapse' bound 

N 

3=1 l<j<£<N 1 3 1 



due to Levy-Leblond ll20l Theorem 2]. Hence the presence of the constant Cq in Theorem ll.2l 
We now turn to the physical regime \/2a < U. Here our lower bound is a byproduct 
of our proof of Theorem 11.21 and we have no reason to believe it is optimal. Together with 
Theorem 11.11 it demonstrates, however, that the model undergoes a sharp transition at a = 

U/y/2. 

Theorem 1.3. For < y/2a < U, 

Vs " J U-V2a 
Last but not least there are the following universal variational upper bounds for En and 

En+m- 

Theorem 1.4. For all N, M, a and U we have 

E N < -aN, 
En+m < En + Em- 

The bound E\ < —a is well known from |fl9l[T0l and it agrees with the result of a formal 
computation of E\ by second order perturbation theory ifTTt . Also, it is consistent with Haga's 
computation of E\ including a 2 -termg3 [16]. The bound En < — aN follows from the esti- 
mates Ei < —a and En < A^^i, the latter of which is a consequence of the second result 
of Theorem 11.41 We remark that En+m < En + Em holds quite generally for translation 
invariant N -particle systems with interactions that go to zero with increasing particle separa- 
tion. In particular it holds for fermions and for distinguishable particles alike. Numerically 
computed upper bounds on E(N)/N, fovN = 2 through N = 32 can be found in the litera- 
ture [5 ], but in the case of fermions they are not refined enough to be consistent with the bound 
En+m < En + Em- 

In this paper we have omitted spin, but the Fermi statistics is taken into account. There 
are only few small modifications necessary for treating fermions with q spin states, such as 
factor of q~ 2 / 3 in front of the Thomas-Fermi kinetic energy, which alters the upper bound in 
Theorem ll.ll by a factor of q 2 / 3 . 

The many-polaron model has also been studied with a confining potential of the form 
Y,eLi W(xi), W(0) = and W > included in the Hamiltonian HH. We could include such 
a potential in our work as well, but, at least in the regime \/2a > U this would not affect the 
leading large iV behaviour of En- 

'There is a sign error in Feynman's quote of Haga's result. 
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2 Upper bounds on E 



N 



In this section we prove Theorem II. II and Theorem 1 1.41 Since En = liniA-^oo En,A we only 
need to deal with the self-adjoint operator H n ,a- Let / G V N = A N L 2 (M. 3 ) n H 1 ^™) 
be normalized and recall that the one-particle density matrix 7 and the density function p 
associated with / are defined by 



-y(x,x'):=N f(x,X2,...,x N )f(x',X2,...,x N )dx 2 '-dx N , (2.1) 

(JV-l) 



p{x) := j(x, x) = JV / \f(x,x 2 ,...,x N )\ dx 2 ---dx N - (2.2) 

In this paper the Fouriertransform /) of the density function p, or of any other function, is 
defined by: 

p(k)= [ e- lk - x p{x)dx, 
that is, without a factor of (27r)~ 3 / 2 . 

Proposition 2.1. Suppose v2a > C/. Then for every one-particle density matrix 7 on L 2 
with < 7 < 1, TV [7] = N, Tr[— A7] < 00, and for p{x) := 7(2;, x), 



£jv < (V2a-U) 2 §TV[-A 7 ]-±/ 

/i 



f -y\ 



x-y\ 
dxdy. 



Proof. This proof is based on the estimate En,a < {f®Vi Hn,aI <S) ? ?) f° r suitable normalized 
f £ T>n and rj E J-. We begin by observing that the expectation value of the interaction oper- 
ator in a state / (g) 77 may be represented in the following two ways: if / and rj are normalized, 
then 



N 



f <g> 77, ®A(xe)f ® 77 

= / |/(a:i,..., xjv)| 2 V"VA.,r,(^)rfa;i---^Ar (2-3) 
= {<n,*A(p)ri) (2-4) 
where VA jr) ( x ) := (*7> &a{x)t]}, p is the density associated with /, and 

*a(p) := / p(x)(»7,*A(a?)i7)daJ 

1 



[^(fc)a(fc) + /3(&)a*(A;)]dA;. 



'|fc|<A c 0| 

Hence if we define if^ A := YleLi[~ + a/^^a,^ 2 ^)] + UVc, then 

(/ ® v, H N ,Af ® 7/> = (/, A /> + (ry, ff ph r?). (2.5) 
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The ground state energy of the iV-body Hamiltonian Hjj A is bounded above by its ground 
state energy in the Hartree-Fock approximation. By Lieb's variational principle, IT23TI and (3j 
Corollary 1], this Hartree-Fock ground state energy is bounded above by 

r.N,A, \ rp I-/ i A . / tt \ i . U f p{x)p(y) -|7(x,y)| 2 

^ jr6 |x — y\ 

for any one-particle density matrix 7 with Trfy] = N and = 7(2;, x). Hence, in view of 
(12.51 ). we conclude that 

En,a < £ H ' F (7, rf) + (r), H ph r]) (2.7) 

for all normalized 77 € T. In order to minimize the right hand side with respect to 77, we use 
that (1231) equals (l2~!4l) . It follows, by Lemma |A!I1 that 



inf [^Tr(V Ai1ll ) + ( V ,H phV )] =~ [ ll y^dk. (2.8) 



7/6^,117/11=1 L ' ^ " J ,/|fc|<A 

By combining (12.61) . (12.71 ). and (12.81 ) and then letting A->oowe arrive at 

E N < \ Tr[-A 7 ] + (U- V2a)\ [ dxdy - % I dxdy (2.9) 

7rs \x-y\ 2 J R6 \x-y\ 

for any one-particle density matrix 7 with Tr(7) = N and p(x) = j(x,x). Here (IA.3I ) and 
(11.21 ) were used also. In the case \/2a = U it is clear from ( 12.91) or from (12.51 ) with r\ being 
the vacuum vector, that En < 0. In the case where /3 := \f2a - U > 0, we choose the 
density matrix 7 on the form 7 = Up^fUt with iTg defined by (Upip)(x) := (3 3 ^ 2 ip(Px). The 
Proposition then follows from UZAUp = (3 2 A and from 7(1, y) = [5 3 ^{l5x : /%) by a simple 
change of variables in the integrals of ( I2.9I ). □ 



The second ingredient for proving Theorem ll.ll is the following lemma. 

Lemma 2.2. Let g € H 2 (R 3 ) with \\g\\ = 1. Then for every p £ L 1 (R 3 ) vv/fn p > ant/ 
J" R3 p(x)dx = N there exists a density matrix 7 such that "y(x, x) = (p * \g\ 2 )(x) and 

Tr[-A 7 ] = I (<ott 2 )1 [ p(x)ldx + N\\Vg\\ 2 . 
Jm.3 



Proof. For the reader's convenience, we recall the proof from 11221 Page 621]. Let M : I 
be defined by M(jp, q) = 1 if \p\ < (6ir) 2j/3 p(q) lj/3 and M(p, q) = otherwise. Then 



(2^)- 3 / M(p,q)dpdq = / p(q)dq = N 
(27r)~ 3 / p 2 M(p,q)dpdq = f (6tt)I / p(q)hq. (2.10) 



We define 7 by 



7= (2^)~ 3 / M(p,q)U pq dpdq 

JR6 

where U pq is the rank one projection given by 



Ilpq^ = 9 Pq / g P q{x)ip{x)dx, g pg (x) = e ipx g(x - q). 
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It follows that j(x, x) = J R3 \g(x — q)\ 2 p{q) dq, and from 

Tr[-An pg ] = \\Vg pq \\ 2 = p 2 + ||V 5 || 2 + 2p ■ (g, -iVg), 
and (12.101 ) we find the asserted expression for Tr[— A7]. □ 

Proposition 12.11 and Lemma 12.21 suggest the definition of a Polaron Thomas -Fermi func- 
tional by 

£ptf(p) := t^(6vt 2 )§ / p(x)§ dx - U E^Mdxdy, (2.11) 
Jr3 J R 6 \x - y\ 

where p G L^M. 3 ) n L 5 / 3 (M 3 ) and p > 0. If p N (x) := iV^iV 1 / 3 ^, then ||pjv||l = N\\p\\i 
and 

7 

£ptf(paO = N 3 £ptf(p)- 
Hence it suffices to consider densities p with f p(x)dx = 1. Let 

£ptf := inf \ £ptf(p) P>0, / p{x)dx = 1 \ 

L JR3 J 

which is finite by Lemma lATTI 
Lemma 2.3. £ptf < 0. 

Proof. Given p G L 1 (M 3 ) n L 5 / 3 (K 3 ) with p > and f pdx = 1, let p B (x) = R- 3 p(R- l x). 
Then J R3 pn(x)dx = 1 for all > and 

£ptf(p R ) = ir 2 y^¥ / P (x)^dx-R- l \f B^Mdxdy. 

Jr3 J R 6 \x - y\ 

This is negative for R large enough. □ 

Proof of Theorem \1J\ Let 5 G L 2 (M 3 ) be given by g(x) = (27r)~ 3 / A e~ x2 / A and set g e {x) = 
e~ 3 / 2 g{x/e), so that \\g £ \\ = 1 for all e > 0. Let = y/2a - U > 0. If (3 = then < 
by Proposition 12.11 Hence it remains to consider the case (3 > 0. Every density function 
p N G L 1 (IR 3 ) with HpaHIi = N is of the form pn{x) = N 2 p(N 1 l 3 x) with ||p||i = 1. From 
Proposition 12. 1 l and Lemma |2~21 combined it follows that 

P- 2 E N <U^ 2 )1( p N (x)Ux-U PN 'j x)pN f {y) dxdy + N\\Vg £ \\ 2 , (2.12) 
7r3 J R 6 \x-y\ 

where p^^ e = pn * \de\ 2 - Suppose 1 < p < 6/5 and let f(k) := \g\ 2 = e _fc2//2 . Then 
PN~ e {k) = p^{k)\g £ \ 2 {k) = p^(k)f(ek) and 



(2.13) 
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By definition of /, by (12. 13b . and by Lemma lATTl 



PN(x)PN(y) dxdy _ [ pnAx)pnM dxdy 



\%-y\ Me \x-y 

\B(k)\ 2 



dk 



1 £ ,-i [ l-|/( g fc)| 2 |/^(fc)| 2 d/; 



\k\* 
']/ 

- 2tt 2 J r3 |/c| 3 ^ 

= iv2 + f e ,-i 2(27r) ,-2^ /■ pMM dxdy . 

7rs f ~ vr 

Combining this estimate with (12.121 ). we see that 

P~ 2 E N < Nh PTF (p)+Ne- 2 \\X7g\\ 2 

+iV 2 +f^- 1 (2 7 r)^ 2 ^ I P {X)P{V) dxdy 

7r6 \x - yf 

for all p G L 1 ^ 3 ) with = 1. If {p n } C L x (]R n ) is a minimizing sequence, SpTF(Pn) —>■ 
Eptf as n — > oo, then 1 1 Pri 1 1 5/3 * s uniformly bounded by (IA.1I ). and hence so is the term 
J p n (x)p n (y)/\x — y\^dxdy for \i < 6/5. Therefore, in the limit n —> oo, we obtain 

0~ 2 E N < mEp TF + \Ne~ 2 + N 2+ 3 e^C^ 

where the constant is finite for p. < 6/5 and ||Vg|| 2 = 1/4 was used. Upon optimizing 
with respect to e we arrive at 

„ 7 9+5 M 

(3 2 E N < N^EpTp + Nz+^Dfj, 

with a new constant D^. This bound with the choice p, = 37/31 < 6/5 proves Theorem ll.il 

□ 

Proof of Theorem \1.4\ We only need to prove that E\ < —a. The bound En < —aN will 
then follow from Ejy-t-M < -EW + Em as pointed out in the introduction. 
Following Nelson [27 ] we introduce 



CO i| fe |<Al(l + f )|jfe| L 7J 

Then 

e lBA H 1A e- lBA = \ (p 2 + 2a* • p + 2p • a + a 2 + (a*) 2 + 2a* a) +F ph -ae A , (2.14) 
where 



k 2 



CO J\k\<A (1 + %- 
1/ 1 

eg J\k\<A \k\ 2 (l + 



e lkx a(k) dk. 



e A : ~ — I E3T ^ 
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From (12.141) we see that, for all normalized / G L 2 

(f ® O, e lB *H lM - lB *f ® n> = </, (-|A)/) - ae A (2.15) 

where SI G J? 7 denotes the vacuum vector. Since inf cr(— A) = it follows from (12.151) that 
£i,A < —aeA, where 

1 /" 1 
lim e\ = — / P5— a& = 1. 

A-oo C 7 M 3 |A;|2(1 + M) 

This concludes the proof of the first bound in Theorem 1 1.41 

A result similar to En+m < En + Em is expressed by Theorem 6 in lfl5l . A copy of 
the proof of that theorem, with small modifications due to the differences of the Hamiltonians, 
also proves the desired bound here. In fact, the main part of the proof of fl5l Theorem 6] is 
Equation (19) and the equation thereafter, which show that the interaction between electrons 
mediated by bosons decreases with increasing particle separation. This part remains valid for 
the coupling function X|/c|<a/( c o|^|) of the Hamiltonian iJjv,A- Other parts of the proof are 
simplified due to the fact the phonon dispersion relation ujlo is constant and hence a local 
operator with respect to the boson position as measured by iVfc. □ 

3 Lower bounds on E N 

In this section we prove Theorems 11.11 and 11.31 The first step is to make sure that phonons 
with large momenta contribute to lower order in N. To this end, for given K, A,8,k > 0, we 
define the operator 

N 



H n ,a,k := -±(l-K)^A, + (l-<5)# ph + [/F c 



N 

/ e -^— \e lk - x *a(k) + e- lk - xe a*(k)} dk. 
t^J\k\<A co\k\ 1 1 



jr[J\k\<A c 

Of course, later on, 8, k £ (0,1) and K < A — > oo. The following result, in the case N = 1, 
is essentially due to Lieb and Thomas ll25l . While a sharp cutoff \k\ < K is used in |[25ll . we 
work with a Gaussian cutoff since we need the Fourier transform of the cutoff to be positive. 

Lemma 3.1. Suppose K,A,a and U are positive, < 5 < 1 and let k := f^^oo* where 
/oo := (x/2-l)/V^. Then 

Proof. For each i G {1, . . . , N}, we introduce three high momenta modes by 

Z f := [ Tf\k)a{k)dk, je {1,2,3}, 



Cq\K\ 
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kj 6 E being the j'-th component of k G M 3 and xa the characteristic function of the set {\k\ < 



A}. For later use we compute the inner product of two functions T- . By straightforward 



computations, 



where 



/ Tf\k)Tf{k)dk = 5 jr ^-U, (3.2) 



Note that 47t/cq = V2/ir and that I^o = liniR^oo Jr = (\/2 — as defined in the 
statement of the lemma. By definition of Hn,a,k, 

N 

H N , A = H NA k + J2(-^ A £ + I k]a)+ SH p* ( 3 - 3 ) 



2 

=i 



if := /" 1 ~ e „f T [e lfc - a!f o(fc) + / t .c.]dfc 
J|fc|<A c oF| 

where we introduced the operators a associated with the ultraviolet part of the electron- 
phonon interaction. The key ingredient of this proof is that 

'^ = t\?^ z f'- z f"] 0.4) 

3=1 

where p^- := —idjd Xgy This identity implies that 

l^lg^l < 2^11^1111(^-^)^11 
3=1 

K 3=1 
3 



K 3=1 

where k > is to be selected, and the estimate 

\(v,(zf?v)\ < \\z®*v\\\\z® v \\ < Uv,(zfzf* + zf*zf) v) 



was used. From (13.21) and Iuk — ^oo it is clear that 



j^iZ^zf + zfzf*) = j^2Zf> Z f + [zf\zf>] 

3=1 3=1 

< ^^ + ^00. (3-6) 



3A" K 
Combining (13.51 ) and (13.61) we arrive at 

TV N 



2 ^ v ; 3kA' p kK ' 
l £=i 



which, by (13.31 ) and the choice « = 8aNI oa /(3K5), proves the lemma. □ 
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Lemma 3.2. Suppose K, A, a, U and n are positive, and < 5 < 1/2. Then 

, , , A A / \/2a \ 2aNK 
H na ,k > - /c) J] A/ - ( - Z7 1 Vc - (3.7) 

Proof. By completing the square in annihilation and creation operators, that is, by using 
Lemma lA72l we see that 

r N r ~K 

(1 - 5)H ph + \e lk ^a{k) + e - lk ^a*(k)} dk. 

C f^J\k\<A |fc| 



N 



fe 

= ?£L_ V / l^e*^-*) dlb - aNK (3 8) 

The integral in (13.81 ) represents the electrostatic energy of two spherically symmetric, non- 
negative charge distributions centered at Xj and X£, respectively, each distribution having total 
charge one, see (1A.3I) . Hence Newton's theorem, [24, Theorem 9.7], implies that 

>M*j-*t) dk < 



\Xj — Xg\ 

Since Cq = 27t 2 a/2, it follows that (13.81) is bounded below by 

y/2a. r aNK 



1-5 ^ (1-8)^?' 
which proves the lemma. □ 

Proof of Theorem \1.2l We shall combine the Lemmas |3. 1 1 and [3721 with suitable choices for 
5 and K. First, suppose that < 5 < 1/2 and that k G (0,1). Since \/2a - U > 0, 
by assumption of Theorem 11.21 the constant multiplying the potential Vc in Lemma 13.21 is 
positive, and hence, after the scaling transformation 

( ^ a ttY 1 h \ 
U 1 — k)x 



1-5 



we may apply (11.61 ) and find that 

> - V inf. ^--A,-y c 



± — ft 

2 



1-k V^2 / V^r 



_ (1=7-^) 7 2aATiv" 
> -Cq- y —N3 - 



1 — K ^/TT 

where Cg is chosen such that (11.61) holds true. 
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We now make the choices 

5 = |jV~i and K = i32/ 00 aiV 1+ i , 

which imply that k, as defined in Lemma [3TT1 obeys k = ^iV" 1 / 9 = 5. Using that (1 — t)^ 1 < 
1 + 2t, for < t < 1 /2, that f7 < \/2a, and /«,/ ^tt = (\/2 - l)/vr < l/(2vr), we find that 

■ffiV A > -C G [\/2a(l + 25)-C/l 2 (l + 2/€)ivI --p-a 2 iV 2+ i 

07T 



-C G (V2a - U) 2 m - (l8C G + — )a 2Ar 

V 37T/ 



> -C G R^2a - C/) 2 + 16a 2 51 (1 + 2k)7vI - — c^ivF! 

> -C G [(V2a - t/) 2 A^ + 18a 2 7vi~5] - — a; 2 ivl-i 

32' 



7_ 1 
3 9 . 



□ 



Proof of Theorem \1.3\ Finally we consider the case, where U — y/2a > 0. In Lemma [3^21 we 
choose S = (U- y/2a)/(2U) and K = SaN 1^/(35), so that k = 1 in Lemma[3j] and 

1-5 2(1-5) 
From Lemma |3~T1 and Lemma 13^21 it hence follows that 

2aNK 3 16Joo 2 2 2C7 3J7 



25 3 " 17 - \/2a U - y/2a ' 

where 7oo/0F = (^/2- l)/vr < l/(2vr). □ 

A Auxiliary Results 

Lemma A.l. Suppose that p G L 1 ^ 3 ) n L 5//3 (R 3 ), p > 0, < p, < 6/5, and let pn{x) = 
N 2 p(N 1 / 3 x). Then 

P -^d xd y < a M l^\\p\\f (A.D 
\x — y^ 3 

P{x)p{v) dxdy = (2vr)-^ f ^dk (A.2) 



\x - # ' ' " " J R3 \k\ 3 

= jl /■ IMP,, (A3) 

F~y| 2vr^ J R3 |fc| 2 

PN(x) PN (y) dxdy = n2+ , f PMM dxdy (A . 4) 



where 



f - yr ' 7ks F - 2/r 



3 / V 5 M/ \5/U / 2 
in (IA.1I ) omg? (IA.2I) , respectively. 
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Inequality (1A.1I) . in the special case p, = 1, implies that £ptf(/o) is bounded below, and 
moreover, that 1 1 >o 1 1 s /3 is bounded unformly for densities p with \\p\\i = 1 and £ptf(#) < 

EpTF + 1- 

Proof of Lemma \AA\ For each R > 0, by Holder's inequality, 



/ 



x-yr j\x- y \<R\ x - vr J\x~ y \>R\ x -y 



By optimizing this bound w.r.to R > we obtain (IA.1I) . Equation (IA.2I ) follows from 
Corollary 5.10]. The factor (2tt)~ 11 stems from the differences in the definition of the Fouri- 
ertransform. Equation (IA.3I ) is the important special case p = 1 from (IA.2I ). and (IA.4I ) is 
straightforward to verify by a change of variables. □ 



Lemma A.2. Suppose f G L Then, for every 5 > 0, 

1 

> — 
- 5 



[<Sa*(ife)a(A;) + f(k)a(k) + /(fc)a*(jfe)] dfe 



2 



a«<i lower bound is attained by the expectation value in the coherent state r] £ T, ||r/|| = 1, 
defined by a(k)r] = —5~ 1 f(k)r]. 

Proof. By completeing the square in creation and annihilation operators 



[8a*(k)a(k) + f{k)a{k) + f(k)a*(k)] dk 



+ — )(a(fc) + — )- — 



2, 

- cZA; 

J R3 L V d / V d J d J 

> -^ll/ll 2 

□ 

Proposition A.3. Suppose that \/2a < U. Then for all N, A > 0, and all f E A JV L 2 (IR 3 ), G 
Twith\\f\\ = h\\ = l, 

40 / 2 \ s 

(f®71,H N ,Af®v) > "CLy ( 

where cl = 1.68 or any of/zer constant for which the Lieb-Oxford inequality holds. 
Proof. As in the proof of Theorem ll.il 

(/ ® T), H NA f <g> 77) 

= (/, (-±A + C7Vb)/) + (77, (F ph + \&MP))V) 

p(x)p(y) 



> A + [/Vc)/)-V2ai / 



dxdy. 
x-y\ 
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Using the Lieb-Thirring ll24l Theorem 2.15] and the Lieb-Oxford inequalities ll22l we find that 

> c LT [ p(x)hx + {U -V2o)h f P ^ P ^ dxdy - c L / p(x)Ux (A.5) 
Jr3 Jr6 \x - y\ J K 3 

where clt = ^g(^f") 2 ^ 3 an d cl = 1.68 or any other constant for which the Lieb-Oxford 
inequality is satisfied. From /?(x) 4//3 = p(x) 5 / 6 p(x) 1 / 2 and the Cauchy-Schwarz inequality, 
for every e > 0, 




(A.6) 



The estimates (1A.5I ) and (IA.6I) with e = 2clt/^cl prove the proposition. □ 
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